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Abstract 



A Newton-Sabatier fixed energy inversion scheme has been used to equate 
inherently non-local p- 12 C potentials at a variety of energies to pion thresh- 
old, with exactly phase equivalent local ones. Those energy dependent local 
potentials then have been recast in the form of non-local Frahn-Lemmer in- 
teractions. 
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I. INTRODUCTION 



Over the years, local nuclear optical potentials have been used predominantly as the 
interaction potentials between colliding nuclei. Frequently, those local interactions were 
specified phenomenologically as Woods-Saxon ( WS) optical potentials; the parameter values 
of which were determined by variation to find best fits to nuclear elastic scattering data. 
The geometries of those WS forms often were taken commensurate with known attributes of 
nuclear densities and/or to be consistent with the character of Hartree fields for the nuclei. 
Most importantly in many cases, quality fits to scattering data were obtained using these 
local potential forms, with parameter values that varied smoothly with energy and mass. In 
reality however, as the optical potentials can be specified in many nucleon scattering theory 
as a folding of pairwise nucleon-nucleon NN interactions with nuclear structure, they must 
be non-local. That is assured by the Pauli principle which implies that the scattering theory 
involves nucleon exchange amplitudes. 

For nucleon-nucleus (NA) elastic scattering, it is now possible to define non-local optical 
potentials in coordinate space and to use them without recourse to localization techniques [Q . 
Solutions of the attendant integro-differential forms of the Schrodinger equations with those 
non-local NA potentials have given good to excellent fits to cross-section and analyzing 
power data for a wide range of target masses and for energies 40 to 800 MeV Those 
non-local NA potentials are formed by folding effective, medium dependent, two-nucleon 
(NN) interactions with the one body density matrices (OBDME) of the target nucleus. 
Very good results for light mass nuclei in particular have been obtained when those OBDME 
were specified from shell model (or similar) calculations involving very large and complete 
shell spaces. 

Nevertheless, for a multitude of uses, such as with application of multi-step reaction 
theories and scattering into the continuum, it is useful to define a local form of optical 
potential for the pA scattering system. Note also that equivalent local potentials can be 
used to account for the coupled channels effects in scattering Q, the resulting form related 
to that found by localization of exchange amplitudes. There are many ways in which a 
local potential can be specified and considered equivalent to a non-local one. Some have 
been reviewed recently JT[|. Herein we consider a form that is to be phase equivalent. The 
phase shifts provided by the g-folding (non-local) potentials used recently M to find very 
good predictions of p- 12 C elastic scattering, have been used as input to finding solutions of 
fixed energy quantum inverse scattering theory. In particular we have found the equivalent 
local potentials by using a modified version of the Newton-Sabatier inversion scheme that 
has been proposed by Lun et al. [[|. This method allows extraction of both central and 
spin-orbit interactions from the input phase shift values. It is not the only fixed energy 
method that can do so. Hooshyar || has used the Sabatier interpolation formulae in a 
finite difference inversion method, and Huber and Leeb J7| have investigated an approach 
to this problem based upon Darboux transformations. Likewise an approximate scheme || 
has been used with some success || to analyze neutron-alpha particle scattering data in 
particular. However, the scheme we adopt is most facile, reducing the process of inversion 
to finding solution of a system of linear-algebraic equations. 

Fixed energy inverse scattering schemes are not the only ways to effect the transition 
from non-local potentials to equivalent local ones and vice versa. Various other approxi- 
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mations and analytic expressions exist [T|]. However, all such other methods rely upon the 
range of the non-locality being small in comparison to the size of most nuclei. If such is also 
the case with regard to the de Broglie wavelength of the projectile, a semi-classical WKB 
approximation may be valid. However, Peierls and Vinh-Mau |T(| note that while the local- 
ization approximations should be well met within the medium of a large nucleus, corrections 
could be important in the nuclear surface region and for light mass nuclei. In any circum- 
stance that the non-locality range exceeds the characteristic length of the system, they note 
that the non-local potential cannot be approximated by their algorithm of localization. A 
desirable feature of using inverse scattering theory to define equivalent local potentials is 
that none of these approximations are invoked. But such must be borne in mind when the 
reverse process, of defining a utilitarian non-local form from the energy variation of a set of 
local potentials, is considered. In doing that reverse process, we chose a mapping used by 
Apagyi et al. |TT| for its simplicity and as the data sets we study exist at disparate energies 
in the range 65 to 250 MeV. 

In this paper the salient features of the modified Newton-Sabatier fixed energy inversion 
scheme are discussed first. Then in Sec. Ill, we present and discuss the results of this non- 
local to phase equivalent local potential scheme. The method to recast non-local forms from 
energy variation of local potentials is then defined in Sec. IV, and results of mapping the 
local optical potentials back to non-local interaction forms of Frahn-Lemmer type [12] are 
given then in Sec. V. Conclusions are drawn in Sec. VI. 



II. THE MODIFIED NEWTON-SABATIER FIXED ENERGY INVERSION 

SCHEME 



The modified Newton-Sabatier method || permits extraction from a set of scattering 
phase shifts, not only of a central local interaction between two colliding quantum objects 
but also a spin-orbit term. The approach relies on analogues of the Regge-Newton equations 
for interactions involving the spin-orbit component of the potential, namely the Sabatier 
transformation equations ||13||. In this method, the non-linear Sabatier interpolation formulae 



essentially are converted into a finite set of linear-algebraic equations which are easily solved. 

With most inverse scattering theories, the assumed equation of motion is the homoge- 
neous Schrodinger equation, the radial terms of which we consider in the form 



d 2 2\iE l{l + 1] 



dr : 



+ 



Rt 







(1) 



V(r) is the hadronic interaction potential which we assume is a sum of central and spin-orbit 



terms, viz. 



V(r) = [V c (r) + V so (r) £ 



(2) 



where cr = 2s. With charged particle scattering, V c (r) also includes the Coulomb potential. 

Since the energy is fixed it is convenient to recast Eq. ([!]) into a dimensionless form, 
using the following notation 



2fiE 



V c (r) 



V„ 



p = kr, k 2 = U c (p) = rc l' ' and U„„ = -^r- 



E 



E 



(3) 
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The dimensionless, decoupled, reduced radial Schrodinger equations subsequently are 



1 
dp 



+ 1 - U c (p) + sU so (p) T 2s\U so (p) 



y- 4 



1 1 



XxU>) = 



(4) 



where A(= £ + |) is the angular momentum variable. 

As two unknowns, V^r) and Ko( r )> are sought, only two of the (2s + l) possible equations 
in Eq. (^) are required. It is convenient to choose the cases for which j = t ± s associated 
with which are two sets of phase shifts denoted by 5f and 8j respectively. This method 
can be used with particles of any spin. For NA scattering, the nucleon intrinsic spin of 
course is |, but we must presume that the spin of the nucleus is not an influential factor 



in scattering, i.e. we equate all such scattering to that from a spin zero target so that the 
quantum number is j = i ± |. This is precisely the case though for p- 12 C scattering that 
we investigate herein. The superscripts '±' now designate the relevant two values of j. 



The Sabatier interpolation formulae relate the regular solutions of Eq. 



Xx(p), to 



the set {ip\(r)}, which are solutions of the Schrodinger equations for a reference potential, 
Uo(r), that are regular at the origin. The link is 



7T 



(p) 



(p) 



(5) 



where 



W x ,(p) 



A 2 -/i 2 

is the Wronskian. The set of scalar functions, F ± , in Eq. (|5D are defined by 

F ± (p) = exp \± [ P tSU s (t)dt 
L Jo 

2 2 

= — aoXo(p)Mp) + — l a txj(p) + b„xt(p)] ^(P) 

71 "P n P fieS 



(6) 



(7) 



and obey the following property: 

F + (p)F-(p) = 1. (8) 

Solution of these equations rely on a complete knowledge of all the phase shifts, 5f, 
where A = £ + 0.5 now includes not only the set of corresponding to the physical set of 
integer I but also all half-integer values. As such, the summations span the set of angular 
momenta f2 : {-, 1, |, 2, |, 3, • • •}. Naturally, analysis of scattering data can provide 
only the physical set of phase shifts, i.e. those corresponding to integer £. Interpolation of 
the (physical) data set is then required to define the required set for inversion. Of course 
if instead one wishes to map against a defined non-local potential, then the phase shifts at 
the intervening values for A can be evaluated. 

Eqs. (|5]) and (0) contain a set of unknown coefficients a^. Eq. (0) also includes functions 
Xo(p) and ipo(p) which are solutions of the relevant Schrodinger equations for £ — — |. In 
addition there are the weights 6a which are defined |14| by 
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7(A + |-^h(A + | + ^) 



7(A + |) 



cosh(7r?7) integer A 
sinh(7rr]) half — integer A 



(9) 



The set of equations, Eqs. (|) and (0), constitute a matrix problem to specify the un- 
known coefficients, a^, and the functions, F ± , when select large radial values pi are chosen 
so that 



*Pf(Pi 
XxiPi 



sin 
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Pi-n(A-~ +aA-77 ln(2pi 



cj" sm 
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77ln(2pi 



F (pi) = h* 1 ] (constants) , 



(10) 



where <ta are the Coulomb phase shifts and are additional unknown constants to be 
determined. Taking two (or more) distinct radial values, (pi = pi,p2,'" > Po), and with 
X\(p) = c \T\(p), Eqs. (H) and (^) can be recast as 



<Ma) = cS=i?(a) + E ^Mpi) fe Pi )Aj - vr±(p;)c, 



(ii) 



They form a set of linear equations in the unknown coefficients which have been grouped as 



Ai 



afc^h T and = c^h^ . 



(12) 



The set S' is limited to X max so that the matrix is finite and we can find a solution by 
using singular value decomposition (SVD); a useful approach since the matrix may tend to 
be ill-conditioned. By so doing we presume that from a characteristic radius, po = kro, 
the interaction is solely Coulombic (or zero for incident neutrons) and that for all A > 
X max (p ), ^a" — <7\ — > 0. With the coefficients af, and defined (/i ± are specified in 
terms of them), multiplication of Eq. (|5]) by F + (p) and upon rearrangement, gives a set of 
(8 * Xmax + 1) linear equations, 

F + (P)] 2 MP) = F+(p) X + x(p) - E -WxM kF + (p)x^(p) - a+F+(p) X ,( P ) 



new 



71 



2p 



MP) = F + (p)xx(p) - E — Wx,(p) b,F + (p) X ;(p) - a;F + (p) X +(p) 



new 



71 



2p 



1 " (F + (p)) 2 UP) = E —WoM (b a + a-)F+{p)xt{p) - ib a + <)F+(p)x M (p) , (13) 



nen' 



71 



that can be solved for the (8 * X max + 1) values of F + (p)x\(p) and F + (p) at each value of 
P < Po desired. Then, as Chadan and Sabatier have shown, the inversion potentials are 
obtained from identities 

u M = ± 2 - dla ™ with F* W r W = i 

p dp 

^(p) = U (p) + if/ so (p) - - p ^- p [G + (p)F-(p) + G-(p)F+(p)} + \[pU so (p)f , (14) 



where 



gHp) = -Y,»HxJ(p) 

71 P iien 



(15) 



Inherent with inversion potentials from all Newton-Sabatier methods is a pole at the 



origin. That pole arises from the s-wave contributions in the summations |T5| but it usually 
influences only small radii properties of the results. In most cases studied, the effects of the 
pole are not evident beyond 0.5 fm typically. 

The Coulomb field poses a problem with inverse scattering theory applications. The 
scattering phase shifts with a Coulomb potential incorporated increase with t. However, a 
transformation of phase shifts allows that problem to be allayed |]IEf| . At and beyond a radius 
r , the nuclear component in the total potential can be ignored. So at r the Schrodinger 
potential is given by V c (r ) = 2r]E/p . The method is to take that value for a new potential 
at all larger radii, and hence that new potential is 



V(r) 



V(r) r < r = p /k 
V c (r ) r>r 



(16) 



For this potential, the internal solutions are then to be matched to Bessel functions (zero 
reference potential solutions). Subtracting the long ranged constant potential V c (r ) from 
the overall one, 



V{r) - V c {r ) 



V c (r) - V c (r ) r <r 
r > r 



(17) 



gives a new potential which when used in the Schrodinger equation gives rise to a new but 
essentially equivalent set of phase shifts. By matching the logarithmic derivatives at r with 
the external, zero potential solutions, one finds the relation 

d ~ d 

— \og[cos5 x Hx(Ppo) + sm5xh(Ppo)} = logfcos 5\i[>x((3po) + sin S x (x(Ppo)} , (18) 

where H\ and I\ are the regular and irregular zero potential solutions respectively and 
(3 = \Jl — V (po) / E cm (= \Jl — %v/Po f° r a Coulomb potential). The new phase shifts are 
then given by 



arctan 



where 



^a(Po) 



' H' x (Pp ) - H x (P Po )D> 
v I x {/3po) - h((3p )D x 

cos5\ip' x (po) + sm8xCx(po) 
cos8 x i/j\(po) + sin5 A CA(po) 



P 



(19) 



(20) 



Thus inversion is made of these new phase shifts with a zero reference potential to obtain 
U(p) from which and on accounting for the energy shift, the actual potential for r < r is 



V{r) = E (3 2 U(r) + 1 - ft 



_P_ 

k(3 



(21) 
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III. RESULTS AND DISCUSSION 



Using medium modified effective NN interactions between a projectile proton and each 
and every bound nucleon in 12 C, microscopic non-local complex and spin dependent optical 
potentials have been generated for a range of energies, 40 to 800 MeV in fact ||. The 
resulting cross sections and analyzing powers so predicted are in good agreement with the 
observed data. Those calculations were made using the program DWBA98 [17]. Therein, at 
each (fixed) energy, the appropriate effective NN interaction is folded with the OBDME of 
the target. For 12 C those OBDME were generated from complete (0 + 2)hu) space shell model 
wave functions 0. That program also solves the integro-differential form of Schrodinger 
equations, and thereby we obtained sets of phase shifts of (non-local) optical potentials to 
use as input to the inversion procedure. 



A. The equivalent local potentials from inversion 

The inversion potentials obtained using the ^-folding model phase shift sets specified for 
proton scattering from 12 C with energies 65, 100, 160, 200, and 250 MeV are displayed in 
Fig. || by the solid, dashed, long dashed, dotted and dot-dashed curves respectively. Results 
have been found for 135 MeV protons as well and they are intermediary to those shown. In 
segments (a) and (b) the real and imaginary parts of the central potentials are displayed 
while segments (c) and (d) respectively contain the real and imaginary components of the 
spin-orbit potentials. The central potential values progressively becoming less refractive and 
more absorptive with increasing energy. The spin-orbit potentials have more variation in 
their structure with energy although there is a general shape for the real and imaginary 
components. 

Clearly while these components of the inversion potentials vary smoothly with energy, 
the central terms especially, their shape is not characteristic of traditional potentials, e.g. 
a Woods-Saxon function and/or its derivative. That is so independent of the inherent pole 
term in the inversion method which dominates the inversion potential values at the origin. 
But, as noted previously JTH], that pole influence is minimal beyond about 0.5 fm. Also at 
each energy, the effect of the pole when the inversion potential was used in calculation of 
scattering observables was small if not negligible. Of primary interest then is the behavior of 
the inversion potentials from approximately 1 fm and our comments pertain to the properties 
from that radius out. 

In Fig. Ill, the real components of the central potentials found by inversion are shown in 
panel (a). At 65 MeV this component varies almost linearly between 1 and 4 fm in radius 
before tapering to zero, and then rather slowly. With increasing energy some structure 
develops in the central real potential for radii 1 < r < 4 fm eventually, at 250 MeV, 
becoming a shoulder shaped well. While values of the real central potential in that region 
(1 < r < 4 fm) decrease in strength to form the shoulder shape, the longer range property 
increases in effect with energy. 

There is a different trend in the imaginary components of the central inversion potentials 
as is evident in segment (b) of Fig. [l|. At all energies the long ranged character of the central 
imaginary interactions (r > 3.5 fm) is the same and the component rapidly vanishes with 
large r, contrasting markedly with the properties of the associated central real terms. With 
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increasing energy the central absorption also increases markedly in the body of the nucleus 
(1 < r < 3.5 fm in this case) with a noticeable structure in the 65 MeV potential result. 
That absorptive character however is almost linear at 250 MeV. 

The real and imaginary components of the spin-orbit potentials found by inversion are 
shown in segments (c) and (d) of Fig. [I]. They are relatively weak, the imaginary parts 
especially so save for the result at 65 MeV. Overall these potentials are very similar; the real 
parts having a weak attractive well (about 2.5 MeV deep at r ~ 1.8 fm) with a shorter ranged 
repulsion. No serious import should be attached as yet to the specifics of the structures 
shown, other than that they are the result of using the particular input set of phase shifts. 
Possible adjustments to the basic NN effective interactions in the original non-local potential 
generation could vary what we are to use as input phase shifts sufficiently to alter the small 
magnitude details of the spin-orbit components shown here. 

B. The Phase Shifts 

Two potentials may be considered equivalent if they give the same scattering phase shifts. 
In principle the inversion potentials will be so. Nevertheless, the inversion potentials must 
be used to get 'inversion' phase shifts for comparison against the original values to check 
that they are indeed equivalent. In Fig. ^| these sets are compared for all energies. The 
real and imaginary parts of the phase shifts values are shown in the left and right panels 
respectively while those for j = £ + 0.5 and for j = £ — 0.5 are given in the top and bottom 
sections respectively. The input phase shift values are depicted by the filled circles while 
the results found from the inversion potentials are indicated by the connecting lines. The 
notation with energy is that used in Fig. |l|. Clearly there is excellent agreement between 
the phase shifts found from the local, inversion potentials and from the ^-folding non-local 
ones. 

There is a relatively smooth and consistent change in the phase shift values as the 
energy increases. The real parts of the sets corresponding to j = £ + 0.5 and j = £ — 0.5 both 
show a steady decrease in value with energy for small partial waves, while the associated 
imaginary parts show a steady increase. The notable exception is the 65 MeV case for 
which the imaginary parts differ from this pattern. In the case of the imaginary parts of the 
j = t + 0.5 phase shifts there is an unusual structure in the 65 MeV data set for the values 
i < 5. 

A closer inspection of the results revealed that the s-wave phase shifts are slightly different 
from the original input values in both the 65 MeV and the 100 MeV cases. This variation, 
however, is less than four percent. The fact that the phase shifts corresponding to £ = 
differ is not unusual since the centrifugal barrier screens all other partial wave solutions from 
any (small) effect of the pole terms inherent in the inversion potentials. 

C. The cross sections and analyzing powers 

Although the excellent agreement between phase shift sets obtained from the local (in- 
version) and from original non-local (full folding optical) potentials for the p- 12 C scattering 
at diverse energies is convincing, another way to demonstrate this equivalence is to compare 
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the associated observables. By so doing any small variation in phase shifts that may exist 
can be emphasized. This is so as the cross section spans several orders of magnitude and 
then small inaccuracies within the phase shifts could be very apparent at the larger scatter- 
ing angles particularly. The analyzing power likewise should be sensitive to small differences 
in the phase shift values since that observable is given by differences between scattering 
probabilities and is normalized by the differential cross-section values. 

The cross sections and analyzing powers at the set of energies chosen (65 to 250 MeV) 
are given in Fig. In the top section of Fig. |3], the cross sections for each energy are 
displayed. Once again the values of the the cross sections obtained from the full microscopic 
NN folding non-local potentials are represented by the circles and squares and the results 
obtained by using the inversion potentials are portrayed by the diverse lines. The coding of 
those lines for each energy is as used in Fig. |1|. Circles and squares, filled and open, have been 
used to display the 'data' to differentiate the set for each energy. These cross section results 
demonstrate that the inversion potentials are very good local equivalents. The cross section 
reproductions span eight orders of magnitude and only for magnitudes less than 10 _2 mb/sr 
are small divergences evident. Such divergence is most evident with the 200 MeV and 250 
MeV cases since those cross sections decrease most rapidly. 

An even finer test of the agreement based upon observables between the local and non- 
local potentials are the results for the analyzing powers. As shown in the bottom section 
of Fig. ^, reproduction of the analyzing power for each energy is very good out to a center 
of mass scattering angle of 60°. Even then only the cases of 65 and 100 MeV have any 
noticeable divergence between the local and non-local potentials results. It is surmised that 
the small variations in the values of low-£ phase shifts are the cause as such variations little 
effect predictions for the cross sections. We conjecture that this behavior is a result of 
the choice we have made for the phase shift value at the (unphysical) angular momentum, 
I = —0.5; a quantity required in the inversion process. 



IV. FRAHN-LEMMER FORMS FROM LOCAL ENERGY DEPENDENT 

POTENTIALS 

The problem is to use local forms of complex potential Vleq(i"; E) = V(r; E) which when 
used in the Schrodinger equation, 

^V 2 0(k, t) + [E- V(t; E)] 0(k, r) = , (22) 
give phase shifts equal to those found by solution of the non-local equations, 

|- V 2 X (+) (k, r) + E X M Qt, r) - J U(r, r'; E) x «(k, r') dr' = , (23) 
where that non-local form is of the Frahn-Lemmer type [1^ . 



If the range of non-locality in U(r, r'; E) is small, then to evaluate the integral term 
in the general form, Eq. |23], of the Schrodinger equation, it is not necessary to know the 
solution function x^ (k, r) at all positions. One only needs to know how x^ (k, r) varies in 
a volume element characterized by a small distance 's' about the point V. In that volume 
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element, x^ + )(k, r) oscillates with a wave number K(r) so that integration over 's' will select 
only those momentum components of any kernel that are in the neighborhood of K(r). This 
justifies expansion of the Fourier transform of a kernel G(s) about the local wave number, 



G(s) 



G(p)e~ tp - S dp 



(24) 



and retention of only the first two terms in the expansion, 



G(p) = G(p 2 ) = G(K 2 ) + (p 2 - K 2 



d 



d(K 2 



-G(K 2 ) + 



(25) 



In a simple manner then a local equivalent potential to the exchange term in Eq. (|23|) can 
be obtained by a Taylor series expansion, 



j U(r, r'; E)x (+) (k, r') dr' » / U(r, r'; E)e i ( r '- r )- v x C+) (k, r) dr' 

= J U (r, r'; E)e i(r ' _r) ' K dr / | x (+) ( k ; r ) 



(26) 



where a local wave number, 

K(r) = ^2/i [E-V(r;E)] , (27) 

has replaced the gradient operator. 

Frahn and Lemmer [12[] assumed that the non-local kernels of the full Schrodinger equa- 
tion, Eq. (0), have a separable form (this is also known in the literature as the Perey-Buck 
prescription), 



U(r,r';E) - F(R)v(p) 



where 



R=-(r + r') 



p = r — r 



(28) 



(29) 



Furthermore, they assume that F(R) = F(R) with F(R) a slowly varying function about 
R = r, and that 



v(p) = v(p) 



1XO 



_3 



exp 



(30) 



where a is the non-locality range. A Taylor expansion about 'r' gives to second order 
F(i?)x«(k,r')« 

F(r) X (+) (k, r) + l -p 2 |^(r)|; + F(r)V 2 + l^[ r F(r)] J X«(k, r) (31) 



and this gives 
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J U(r,r';E)x {+ \k, r') dr' 



voF(r) + g«2~ [rF(r)] + \^/{r)^ + ^(4 (+) (k, r) 



where t> n are the moments of the non-locality. 

47T 



v(p)p {2+n) dp 



n (2n - 1) Jo 

Under this approximation, the non-local Schrodinger equation reduces to 



4 dr dr 



-E + F(r) + a 2 -^ [rF(r)}\ X (+) (k, r) = 
r dr z J 



which maps to an equivalent local and energy dependent form, 

h 2 



-L-V 2 + V{r,E)-E 

— Lb 



¥>(k,r) =0 



upon using a point transformation defined by 

Vp(k,r)=T(r) X (+) (k,r). 
Multiplication of Eq. (Rl) by T(r) gives 



T(r) 

with the functions, 



j-^v 2 - 


\Ji(r)~ 


d 


~Ur) + E~ 


! 


[X(r)_ 


dr 


[ X(r) _ 





|T- 1 (r) V 9(k,r)=0 



J (r) = F(r) + 



a 2 d 2 



16r <ir 2 



[rF(r)] 



X(r) 



2^ 



F(r) . 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



The first derivative term is eliminated by choosing T(r) = JX(r) as then Eq. ( |37D becomes 

T(r) — E' 



V 2 

2/i 



X(r) 



so identifying 



V(r,E) = E + 



<p(k,r) = 



T(r) - £ 



(39) 



X(r) 



(40) 



11 



when 



Mr) 



3 1 



2X(r 



a 



2 d 



8r dr 



a 



2 d 2 



16r <ir 2 



Ffrl 



3/iO" 4 



32/T 



dr 



(41) 



Following Apagyi ei a/. ||11|| , by considering local equivalent potentials at two energies 
Ei and £7 2 , it follows by using Eq. (HO) that 



T(r) 



and that 



F(r) 



h z 2 



E 2 V(r, E 1 )-E 2 V(r,E 2 ) 
[V(r,E 1 )-V(r,E 2 ) + E 2 -E 1 ] 

V(r,E 1 )-V(r,E 2 ) 



H a 2 [V(r, E^ - V(r, E 2 ) + E 2 - E±] 



(42) 



(43) 



Then from Eqs. fl38|) , one can find Jo(r), Ji(r), and -X"(r), and so have a complete specification 
of the non-local interaction properties. To the extent that the energy dependent Frahn- 
Lemmer form describes p- 12 C scattering the functions, U(r), F(r) and T(r) should not be 
energy dependent. As will be seen, that is not completely the case with the system we have 
studied. 



V. THE FRAHN-LEMMER POTENTIALS FOR P - 12 C 

The central parts of the inversion potentials that were specified starting with phase shift 
sets from the g folding optical potential calculations of p- 12 C scattering at 100, 135, 160, and 
200 MeV have been used in this study. With the non- locality range, a, taken first as 0.7 and 
subsequently as 1.0 fm, the pairs of inversion potentials with energies 100 and 135 MeV, with 
135 and 160 MeV, and with 160 and 200 MeV, have been used to find the functions U(r), 
F(r), and T(r) that characterize the non-local Frahn-Lemmer form of the p- 12 C optical 
potential. Results found using those pairs are identified by the notation 100-135, 135-160, 
and 160-200, and are portrayed in the next three figures by the solid, long dash, and short 
dash curves respectively. The real and imaginary parts of the various functions are given in 
the top and bottom sections of these figures with the results found using a = 0.7 and 1.0 
fm presented in the left and right side panels respectively. 

The results for the local attribute U(r) of the Frahn-Lemmer representations for the 
non-local optical potential are given in Fig. f|. In general those components are similar for 
all three energy pairs, and more so for the 135-160 and 160-200 cases. The variations in 
these results are but a few MeV in size. However there is a noticeable change in the degree 
of structure with marked oscillations in the results found with the smaller (0.7) non- locality 
range. Still the results with both non-locality range values do exhibit a residual energy 
dependence, and such is presumed not to be the case with the Frahn-Lemmer prescription. 
The modulating functions, F(r), of the actual non-local term in the Frahn-Lemmer form 
are displayed in Fig [5]. The real parts of this function are quite similar in the 135-160 and 
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160-200 cases although the overall strength of the real part of F(r) decreases with the energy 
as it does with increase in the non- locality range. The 100-135 results (for the real part of 
F(r)) are larger than the others and varies from those in structure. The imaginary parts 
of F(r) vary noticeably with the marked structure of the 100-135 results diminishing with 
energy. But the size change is not linear. Such energy dependence is also at odds with 
the Frahn-Lemmer requirement of an energy independent non-local potential. The local 
elements T(r) for these cases are shown in Fig. They strongly reflect the properties of 
the relevant F(r) albeit that the structures are enhanced. 

In the next three figures the 100-135 and the 160-200 inversion potential pairs are shown 
again but now to compare more directly the effects of different choices of the non-locality 
range, a. Results found with a = 0.7, 1.0, and 1.4 fm now are displayed by the solid, 
long dash, and short dash curves respectively. Again the real and imaginary parts of the 
characteristic functions U(r), F(r), and T(r) are given in the top and bottom sectors of the 
diagrams. The U(r) components are displayed in Fig. [7|. The real parts for both energy 
pairs are similar with a decrease in the structure of the results being evident as the non- 
locality range increases. Indeed the results with a = 1.0 and 1.4 fm are very similar. The 
imaginary parts of U{r) have the form of an attractive well. Again the structure observed 
with the shortest range washes out when a is increased. Now also the 160-200 MeV result 
changes strength noticeably with increase of the non-locality range. The results for F(r) 
and T(r) are given in Figs. || and || respectively. The T(r) variations reflect those of F(r) as 
before. In these cases changing the non-locality range has a dramatic effect, largely upon the 
magnitudes. In part though, that might be considered just an off-set to the normalization 
which depends on a. But there are also changes in the structures. The disparity between the 
100-135 and 160-200 results emphasizes again that there is a residual energy dependence one 
must consider if the set of fixed energy inversion potentials are recast as a Frahn-Lemmer 
type of non-local interaction. 

VI. CONCLUSIONS 

The fixed energy inverse scattering method of Lun et al. || has been used to specify 
local potentials from sets of phase shifts given by solutions of the Schrodinger equations with 
non-local optical potentials for proton- 12 C elastic scattering for a range of proton energies. 
Those (non-local potential) phase shift sets give very good fits to observed cross-section 
and analyzing power data M. So also then do the local (inversion) potentials as we have 
shown them to be both phase shift and observable equivalent (to the non-local potential 
expectations) to a high degree. As a non-local to local potential conversion scheme, the 
inverse scattering theory method has proved to be very effective. Irrespective of the pole 
term inherent in the inverse scattering theory method, the local inversion potentials do not 
resemble the simple functional forms, e.g. Woods-Saxon potentials, that are commonly used 
in phenomenological (numerical inversion) analyses of such scattering data. 

We then considered a reverse mapping to see if and how the energy dependence of the 
inversion potentials might reflect a non-locality of simpler functional kind, and of the Frahn- 
Lemmer form in particular. The results are indicative of a characteristic forms for the diverse 
components of that simple non-local form of the interaction, but the detailed properties can 
vary significantly with the choice of the non-locality range, and there is an energy dependence 
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residual in the Frahn-Lemmer functions. Thus we contend that the energy dependence of 
local potentials for p- 12 C scattering (and by implication for other targets) whether those 
potentials are found by inversion or by phenomenology, is not solely a reflection of the true 
non-locality in the interaction between the nuclei. The non-locality itself is also energy 
dependent. 
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FIG. 1. Inversion potentials from p- 12 C scattering analyses portraying (a) the central real, (b) 
the central imaginary, (c) the spin-orbit real, and (d) the spin-orbit imaginary components. The 
separate energy results are identified in the text. 
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FIG. 2. The phase shifts in radians (real parts on the left, imaginary parts on the right) 
obtained using the inversion potentials given in Fig. || compared with the values used as input to 
the inversion procedure (filled dots). The lines connecting the physical values (integer I) identify 
the energies with the same scheme as used in Fig. |]. 
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FIG. 3. The differential cross sections (top) and analyzing powers (bottom) obtained from use 
of the inversion potentials of Fig. |l| compared with the values associated with the phase shifts used 
as input to the inversion procedure. The lines indicate the disparate energy values as specified in 
the text as do diverse symbols for the 'data'. 
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FIG. 4. The local components, U(r) of Eq. (|23|), obtained with non-locality ranges, a, of 0.7 and 
1.0 fm (left and right panels respectively) and as deduced from the 100-135 MeV (solid curves), 
the 135-160 MeV (long dash curves), and the 160-200 MeV (short dash curves) pairs of (local) 
inversion potentials. 
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FIG. 7. The local terms U(r) found using the 100-135 MeV (left) and the 160-200 MeV (right) 
pairs of inversion potentials but for Frahn-Lemmer non-locality ranges of 0.7, 1.0, and 1.4 fm. The 
results are portrayed by the solid, long dash, and short dash curves respectively. 
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